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Aim of this work is to calculate explicitly the result of the experiment of superposition of a mirror 
in the Michelson photon cavities interferometric device proposed by Marshall, Simon, Penrose and 
Bownmeester, as expected within a recently proposed model of non-unitary self-gravity inducing 
localization. As for other proposals of modifications of Quantum Mechanics in a non-unitary sense, 
aimed to account for both unitary evolution and irreversible collapse, like in the famous Ghirardi- 
Rimini-Weber and Pearle's models, it turns out that, for the experimental parameters proposed, 
no effect is detectable at all. It is pointed out that the enhancing properties of matter granularity 
does not substantially change this conclusion. Parameters have also been exploratively varied in a 
certain range beyond the proposed values. It is shown that within 'sensible' parameters, that are 
not yet attainable within current technology, the model exhibits a peculiar signature with respect 
to other collapse models as far as parameters space is explored. Besides, the calculation offers a way 
to see non-unitary gravity at work in a quasi-realistic setting. 
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I. INTRODUCTION 

In the past few decades several proposals of modifica- 
■tion of Quantum Mechanics (QM) have appeared in the 
literature, aimed at unifying its internal fundamental di- 
chotomy between unitary deterministic quantum dynam- 
ics and non-linear irreversible state collapse following a 
measurement process On the other hand big efforts 
have been devoted towards an attempt to reconcile Ein- 
stein gravity with quantum theory. In this context, some 
approaches have focused on the possible role of gravity 
in state function collapse as a result of the incompatibil- 
ity of general relativity and the unitary time evolution 
of QM It has been shown, in fact, that the exis- 

tence of linear superpositions of states with macroscopic 
mass-distribution differences would entail a breakdown of 
classical space-time making the traditional quantum dy- 
namics somehow troubling Q . As distinct from Penrose 
proposal, some other detailed collapse models have been 
proposed, which are based on a spontaneous stochas- 
tic state vector reduction: the Ghirardi-Rimini- Weber 
(GRW) model Q, the quantum mechanics with univer- 
sal position location model (QMUPL) of wave function 
collapse [H and the continuous spontaneous localization 
(CSL) model Recently the mechanism of spontaneous 
symmetry breakdown of time translation symmetry has 
been invoked as well in order to give rise to the quantum 
state reduction [f| . 

In a different proposal, De Filippo introduced a 
nonunitary model of Newtonian gravity (NNG, from now 
on), which can be seen as the non-relativistic limit of a 
classically stable version of higher derivative gravity (see 
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e.g. Ref.Q, references therein). This model presents 
several appealing features to become a natural candidate 
for an effective low-energy model of gravity. For exam- 
ple, while reproducing at a macroscopic level the ordi- 
nary Newtonian interaction, it presents a mass threshold 
for gravitational localization, which for ordinary matter 
densities is around fO 11 proton masses [TcJ. The model 
can be seen as a realistic version of the nonunitary toy 
models [IT| - [l3j inspired by the emergence of the infor- 
mation loss paradox |14l - tl6| from black hole physics. On 
the other hand the violation of unitarity when match- 
ing quantum mechanics and gravity was argued also out- 
side black hole physics, on general consistency grounds 
0, [f7| • The model affords a mechanism for the evolution 
of macroscopic coherent superpositions of states into en- 
sembles of pure states, each one of them corresponding - 
within a future consistent general covariant theory - to 
an unambiguous space-time. Its features include its abil- 
ity to produce an evolution of the density matrix com- 
patible with the expectations leading to the phenomeno- 
logical spontaneous localization models, as it was argued 
that they should be both nonlinear and nonunitary 
While sharing with the other proposals the non-linear 
non-unitary character, at variance with them, however, it 
does not present obstructions consistent with its special- 
relativistic extension [l8[ . Another success of the model 
is the emergence of a unified picture for ordinary and 
black hole entropy as entanglement entropy with hidden 
degrees of freedo m |9| , in agreement with Bekenstein- 
Hawking entropy [f9( and Hawking evaporation temper- 
ature; that arises from the smoothed singularity of the 
black hole introduced by the model and paves the way 
for the quantum foundations of the second law of ther- 
modynamics. 

It is important to realize that the subject of a fun- 
damental non-unitarity, and the various detailed mech- 
anisms proposed 0-(3, [DI I2fjl - l22j |. is not just a matter 
of philosophy, but could be, in principle, experimentally 
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proved or disproved. Current technological progresses 
which have being achieved in isolating, manipulating and 
controlling an higher and higher number of degrees of 
freedom indicate a not far possibility of detecting fun- 
damental decoherence, which would manifest in a clean 
way only once the system has been sufficient ly p rotected 
against the sources of environmental noise |23| . Indeed 
an experiment designed to detect fundamental deviations 
from unitary quantum evolution would be of consider- 
able importance. Some technologies and devices have, 
at present, been recognized to be particularly suitable 
to create quantum state superpositions which are macro- 
scopically distinct (24|. Among them, there are diffrac- 
tion of complex molecules up to 2 x 10 3 proton masses 
[25l | , current cat states in SQUID devices (2(| and super- 
positions of atomic matter waves in Bose Einstein con- 
densates |27l |. 

Recent progress in optomechanical systems may soon 
allow one to make superpositions of even larger objects, 
such as micro-sized mirrors or cantilevers [281 ] . and to 
test quantum phenomena at larger scales. In this con- 
text an appealing quite recent proposal for the practical 
realization of the Penrose "gedanken experiment" @ con- 
siders the relatively small CM-displacement of a lump of 
10 proton masses in a interferometric device in which 
two high-finesse optical cavities are inserted into its arms 
[29| . The cavity in arm A has a very small end mirror 
mounted on a micro-mechanical oscillator (cantilever), 
which suffers the radiation pressure of the photon con- 
fined inside it and as a consequence can be excited into 
a distinguishable quantum state. A single photon inci- 
dent on a 50-50 beam splitter will realize a superposition 
of being in either of the two arms; then, the coupling 
between the photon and the cantilever will lead to an 
entangled state putting the cantilever into a superposi- 
tion of distinct positions. After a full mechanical period 
of the cantilever, it recovers its original position; if the 
photon leaks out of the cavity at this stage, a revival 
of the interference (visibility) is observed, provided that 
the quantum superposition state of the system survives 
at the intermediate times. Conversely, if the state of the 
system collapses due to some decoherence mechanism, 
visibility will not revive. Summarizing, a measurement 
of the magnitude of the revival of visibility gives a mea- 
surement of decoherence occurred in the time interval 
under consideration. 

Our work is devoted to calculate explicitly the output 
of this experiment [29( according to nonunitary Newto- 
nian gravity model [J [1(3, H3, El ■ We would point out 
that detailed calculations for the expectations of some 
other collapse models, gravitational or not, have already 
been done, demonstra ting a far reaching possibility to 
confirm the theory [32h34| . 

The plan of the paper is as follows. Section II contains 
a brief description of the basic NNG model. In Section 
III we give a qualitative discussion to argue, in the rel- 
evant parameter space of the experiment by Marshall et 
al. [20], the subspaces where gravitational effect could, 
in principle, be visible with both homogeneous and gran- 
ular assumption on mirror mass distribution. Section IV 



is devoted to the application of NNG model to the exper- 
iment, and its general solution. Then, in Section V, we 
use the Wigner function to monitor the mirror's state, 
and verify its behavior in time, after a measurement of 
photons' state. Finally, in Sec. VI, we draw conclusions 
and outline perspectives of this work. Calculational de- 
tails are devoted to Appendices. 



II. NON-UNITARY NEWTONIAN GRAVITY 
MODEL 

The aim of this Section is to briefly recall the key fea- 
tures of the NNG model which we will use for calcula- 
tions. On the basis of a number of considerations (among 
which consistency with basic formal relations of QM, ap- 
proximate energy conservation, classical and quantum 
behavior of matter, requirement that non-unitary terms 
have a gravitational origin, etc), it is possible to isolate 
a two-parameter class of non-unitary gravity models, as 
discussed in detail in Refs. [§, [ltj Ho, H3- We will com- 
ment later in the Section on these parameters, while we 
give here a concise definition of the model in its simplest 
form, which will allow us to carry out our calculations. 

Let H[ip\ V] be the non-relativistic Hamiltonian of a 
finite number of particle species, like electrons, nuclei, 
ions, atoms and/or molecules, where tp T, tp denote the 
whole set ipj(x), tpj (x) of creation-annihilation operators, 
i.e. one couple per particle species and spin component. 
H[ijjt,ip] includes the usual electromagnetic interactions 
accounted for in atomic, molecular and condensed-matter 
physics. To incorporate that part of gravitational in- 
teractions responsible for non-unitarity, one has to in- 
troduce complementary creation-annihilation operators 
ipj(x),ipj(x) and the overall (meta-)Hamiltonian: 

H tot = H[^,TP]+H[ft,4,} + 

„v- /* 1 j ^](x)Vj(x)^I(y)^(y) 

- G } nijirik / dxdy j j , 

(hi) 

acting on the product ® of Fock spaces of the 

ip and ip operators, where m, is the mass of the i-th 
particle species and G is the gravitational constant. The 
%p operators obey the same statistics as the corresponding 
operators ip, while [ip,tl>]- — [ip,tl>']- = 0. 

The meta-particle state space S is the subspace of 
Fij) <8> F^, including the meta-states obtained from the 
vacuum 1 10)) = |0)^, ® |0)^ by applying operators built in 
terms of the products ipj (x)ipj (y) and symmetrical with 
respect to the interchange ip* -f-> ift; as a consequence 
they have the same number of ip (physical) and ip (hid- 
den) meta-particles of each species. Since constrained 
meta-states cannot distinguish between physical and hid- 
den operators, the observable algebra is identified with 
the physical operator algebra. In view of this, expecta- 
tion values can be evaluated by preliminarily tracing out 
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the ip operators. In particular, the most general meta- 
state corresponding to one particle states is represented 

by 

||/» = [dx I dy/(x,yW(x)$( y ) 1 1 0> > . (IL2) 



/(x,y) = /(y,x) 

This is a consistent definition since H to t generates a 
group of (unitary) endomorphisms of S. 

If we prepare a pure n-particle state, represented in 
the original setting, excluding gravitational interactions, 

by 



|.g) = j d"x. 9 (x 1 ,x 2 ,...,x n )4(x 1 )v4(x 2 )...^(x„)|0) 
its representation in 5* is given by the metastate 



(Tx (Pylg (x x , x n ) g (y x , y„) x 

x y4( Xl )...4>„) ^t i ( yi )...^ (y j|| 0)) 



A comment is in order on the possible extensions of 
the model outlined in Refs. said in the 

beginning of the Section, the phenomenological general 
model depends on the two parameters (N,e). The first 
refers to the number N of copies in interaction, which on 
thermodynamical grounds can be inferred to be 2 (as in 
the model presented above). It is however interesting to 
note that the limit N — > oo (with e = 1) reproduces the 
famous non-linear Newton-Schrodinger equation, some- 
times considered in the literature as a possible candidate 
equation for the inclusion of the self-gravity in QM, rele- 
vant to the quantum-classical transition [Io| . The sec- 
ond, e, modulates the degree of nununitarity encoded in 
the gravitational interaction. The above model definition 
corresponds to e = 1 , for which all Newtonian interaction 
is of nonunitary type. This choice has been made in order 
to maximize the effect of nonunitarity, while calcula ting 
the prediction on the experiment by Marshall et al. [29| 
in the best model-setting which gives the largest possible 
deviations from unitarity. 



III. A SEMI-QUANTITATIVE ARGUMENT 
FOR NNG EFFECTS IN MIRROR EXPERIMENT 

Before considering the detailed application to the mir- 
ror experiment of Ref. [2j|, in this Section we give a 
semi-quantitative argument for a gross identification of 
NNG effects. 

When considering self-interaction gravitational energy, 
the threshold mass of localization, M tr ~ 10 11 proton 
masses (= 1.672 x 10~ 16 Kg), can be identified in the 
following way. If M < M tr the metasystem behaves like 
an hydrogen-like system, while in the case M > M tr 
the hidden mass-copy is quite well superposed to the 



physical one. As a consequence, the interaction potential 
can be approximated, within the lowest energy state of 
the meta-system, by the harmonic oscillator ground state 
with gaussian wave function width 



1/2 
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where p s u — 5x 10 3 Kg/m 3 is silicon density. For nonuni- 
tary gravity to be effective in localizing the mirror, this 
length scale must be at least comparable with the wave 
packets separation: 



Ax 



2Mw r , 



(IILl) 



where Lo m is the mirror's frequency and k is the optome- 
chanical coupling constant; then the condition Acc > A^ 
amounts to 



1 (u m \ : 

PsU G \ K 2 ) 



< 1; 



for the experimental parameters, i.e. 
500Hz, p sii ~ 5 x 10 3 Kg/m 3 and 



(III.2) 

~ 2tt x 
1, we get 



Ax ~ 5.79 x 1CT 14 to, x ~ 10". As pointed out in Ref 
[35l | , an enhancement in the possibility of observing grav- 
itational decoherence effects is provided by taking into 
account the real distribution of mass inside a crystal, 
which is very concentrated within nuclei (see Appendix 
A. 2). In that case, one should consider instead of p s u a 
matter density p nuc ~ 10 p s u, given by a silicon nucleus 
mass divided by its effective radius, which can be esti- 
mated as the typical spread of the wavefunction inside a 
crystal. This leads to x ~ 10 9 , which is still much greater 
than unity. 

It should be stressed that the choice between homoge- 
neous or granular matter distribution is not arbitrary, but 
is dictated by the experimental situation. As a matter of 
fact, when the relative displacement of meta-masses Ax 
is of the order of the nucleus effective radius, it seems 
appropriate to take into account granularity; when Ace 
is made much greater than interatomic separation then 
homogeneity assumption appears to be the most suit- 
able one; finally, for Ace of the order of interatomic sep- 
aration, if imperfections like dislocations are present in 
the sample (as it usually happens, even when very ac- 
curate preparation methods are used) , then meta-masses 
are likely to 'feel' an effective homogeneous masses poten- 
tial (this is because in the presence of a sufficient number 
of dislocations, as two meta-nuclei get nearer and nearer 
in one place, two nuclei in another place in the crystal 
can equally well go farther and farther from each other) ; 
otherwise, for a really perfect crystal granularity should 
come again into play. 

On the other hand, (fundamental) decoherence rate 
must be at least comparable with (or lower than) a period 
of natural oscillation of the mirror, which, in its turn, 
must be comparable with (or lower than) environmental 
decoherence rate for the experiment to be feasible. 
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Then the following chain of relations must be satisfied: 



It ac J 

h 3w m 



> 



Id, 



(III.3) 



where Eq is the gravitational interaction energy of the 
meta-masses (see Appendix A), is the environmental 
decoherence rate of the mirror [29, 32] . With the parame- 
ters of the experiment the value w m = aj^ p ~ 2ir x 50QH z 
has been proposed. 

An exploration of parameter space within the exact 
solution has confirmed that first inequality (|III.2[) and 
approximate equality in Eq. (|III.3[) must hold in order 
to see some relevant deviation from unitarity. 

It should be stressed that, in spite of the improvement 
of mass size in the mirror experiment with respect to 
double-slit type experiments, the degree of macroscopic- 
ity of superposition is controlled not only by mass but 
also by space separation of the superposed wave packets. 



A comment is in order on the apparent independence 
of the above result on the mass. The main point is that 
we have chosen to measure CM state separation in units 
of coherent states' size, which means that a larger mass 
is associated with a smaller unit of length and then, fix- 
ing k, to a narrower peak separation. One could also 
choose to fix the absolute separation Ax, and express 

k = k (M, L0 n 



IV. 



A. 



2Mlu„ 



in the above formulas. 



APPLICATION OF NNG MODEL TO THE 
MIRROR EXPERIMENT 



In this Section we carry out interference visibility cal- 
culations within the mirror experiment [29j j in the frame- 
work of the NNG model [1, [lO, H3, ISH . 
To this aim, let us start by defining the gravity-free 
Hamiltonian 



J 



Hf ree [b, 6 f , Af A , M B ; u> m ] = hjj Ph (Na+Nb) + hu m tfb - hgM A (b + 6 f ) 



r 



where g = mo m , Ma,b are the number operators for the 
photon in the interferometer arms A and B respectively, 
while b and 6^ are the phonon destruction and creation 



operators associated with the motion of the mirror's CM. 
In this way our (meta-)Hamiltonian can be written as 



H,, 



H free [b, b\N A ;Lo* m ] + H free b, b\N A ; < - K G {b + fo+) (6 + 6+) 



r 



with uj* n = ujm^jl + 2j£f- (see Appendix A for the cal- 
culation of the gravitational interaction strength Kq in 
both homogeneous and granular case). In practice, the 
(relevant degrees of freedom) meta-system is formed by 



two gravitationally coupled harmonic oscillators and two 
by two photonic modes, each couple of modes interacting 
with its own mirror. Then, Schrddinger state at time t 
is given by: 



J 



|*(t))) = ||*(t))) Soh 



d 2 pd 2 f3 K M ^ A 08, p- i)|V(«)>/3 ® m))p , 

I 



(IV.l) 



where 



Itf (*)>„ = 
1 



lUpht (\0a1b) ® |/?c) + fW) IUOb) ® |A)) • 



with 

f(/3) — e iK2 ("™ t ^ sinLJ .*r 1 *) e iK3 '[' 3 ( 1_e 



-iw* t 



Pe 



nil — e 



— iuj* t 



Here Q(x) denotes the imaginary part of x. Computa- 
tional details on kernel K are devoted to Appendix B. 

Since the only experimentally accessible quantity is the 
visibility, defined as twice the off-diagonal term (in the 
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absolute value) of the physical photon's state pab, we 
are going to calculate this quantity as 



Vis (t) = 2 



7Y _ Tr~ 



It can be shown that visibility has a one-to-one corre- 
spondence with von Neumann entropy, which represents 
a good measure of entanglement for a pure bipartite state 



where Ry is defined as 



R ( v ] = (U0b||*»«*||IU1j 



We write the kernel as: 



K"*"*(j3,0,t) = A(i)e 



2h 2 



4a I r 1 * ' r+ +a3=f -1* ' r + +F T 2 - T - 



x e 2 
x e 



x e 



-+/3*a+/3*a 



-7,r_ 



where the functions Ti = Tiiuj^J) are defined in Ap- (fc' = Kq/2Ti 2 ): 
pendix B, Eq. (|B.3|) . Then we proceed to write the 
products K MA ^ A (P,p,t)K M ^*(0',P',t) in the form 



^MaMa t t \^M' A M' A * n\ e -k'[F p l3' 2 +2F rp l3*} -k' [F p 0* 2 +2F m p* j3* +2F Tm /3*] £ -k' [f; /3' 2 +2F* p f3'} -k' [f; /3' 2 +2F^/3' 0' +2F* m /3'} x 
W Wv » ' v » * , '" , ' 



Let's calculate the traces: 
y r _ y, — d(°) 

- 1 ' m,m J 7 p/ l Jt V r 



4tt 4 



9(/3,/3') 



(cj* i— sin cu* t) », 



+ K 11 (f3,j3)K 01 * (ft ,/?') e^I^i-e-^-')] g-^^'ti-e-™™')] e ftft'*_i|ft| 2 _i|^'| 2 



<5(/3') 



Hi(P,P') 



Visibility is then given by two contributions: 



Vis (f) 



7> „ rp r>{a) 



2|(/) + (I/)|. (IV.2) 



Here we do not write the explicit integrals for (/) and 
(//), which can be found in Appendix C; it is worth not- 



ing that, if we discard their photons number dependence, 
they are formally similar. In the free case of no NNG in- 
teraction, i.e. by putting k' — 0, and for a = 0, we 
get 



G 



(J) = (II) = I|A(i)|V K2( " mt - sin " mt) e- K2 ( 1 - cos "'" t ) (for k' = 0, uj* m = uj m and |A(t)| 2 = l) 



as found in 12 i 
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FIG. 1: Interference visibility as a function of time in the 
homogeneous case plotted for k = 1, mass = 5 x 
for the mirror size L = 10~ J m and several values of the fre- 
quency: uj m = x 10~ 5 (dot-dashed), u m = X 10~ 6 
dashed, uj m = 5oj^ p x 10" 7 (dotted) and cj m = x 10~ 7 
(continuous black line). Visibility and time are expressed in 
dimensionless units. 



The behavior of visibility for k = 1 and four different 
values of uj m is depicted in Figure [T] The first case shows 
no difference with the free case, and the remaining three, 
in decreasing order, progressively give a more and more 
reduction of revival effect at the end of the cycle. A direct 
comparison can be made with the predictions of other 
phenomenological collapse models, like GRW, QMUPL 
and CLS @,[3|. Let us focus, for example, on the GRW 
model ■ In such a case the explicit formula for visibility 



Vis (t) = Vis (t) 



e 2 ft. 



-B '10 



where Visq is the visibility in the free case, i.e. in the 
absence of any mechanism of decoherence (of fundamen- 
tal or environmental nature), r/o ~ 0.5 x 10 _2 s _1 m -2 
and Unuc ~ 1-67 x 10~ 27 kg is nucleon mass. Notice 
that also in this case no explicit dependence on mass 
appears. It turns out that in the "worst" NNG case 
(cj m = uj°* p x 10~ 7 ), in which visibility goes practically 
to zero, GRW predicts a visibility lowering at the end of 
the cycle of about 0.5% with respect to the free case. It is 
interesting to note that for smaller Aa; an enhancement 
of the effect's observability comes from the consideration 
of granularity. As discussed in the previous Section, in 
fact, for Aa; ~ 10 -12 or less matter granularity may come 
into play; as a consequence, we obtain that already for 



uj m = lj^ p x 10~ 3 visibility behaves in a similar way to 
the homogeneous case with oj m = cjj^ p x 10~ 6 (dashed 
curve of Figure [IJ . 



MONITORING THE MIRROR'S STATE: 
THE WIGNER FUNCTION 



In the previous Section we have considered a mea- 
sure of photons' interference, and shown that with a 
proper choice of parameters one gets a reduced revival 
effect at the end of one oscillation period. Here we want 
to elucidate about the mirror's state soon after a pho- 
ton measurement process. For monitoring mirror's state, 
and to get a physical insight of what's going on, we use 
the Wigner function. As it is well known, this quasi- 
distribution have both positive and negative parts, the 
latter being a signature of quantum coherence survival. It 
is expected that the action of NNG-induced decoherence 
would reduce, after some time, the interference patterns 
in the Wigner distribution. 

We calculate the Wigner function starting from the 
expression: 



W(x,p;t) 



with 



1 



2Tm 2 



Prn (t) e 



(v.i) 



ip 



2h ' 



where p m is the reduced density matrix of the mirror af- 
ter a photon detection [3a]. Following [32|], we take this 
measurement as the process projecting (physical) pho- 
tons' state onto the state 

\^ = -^={\0 A ,l B )+e* 6 \l Al B )), 

where 8 is a phase constant (it can be shown that the cor- 
responding results are quasi-independent of 9 and then 
we put 8 = 0). Correspondingly, mirror's density matrix 
iB(\\<p)) = \<p)®\!p)) 

p m (t) = Tr ph Ph fri [(Hp)) (M| ®1) ||* (t)» «* (t)\\] . 
After some calculations, we get the Wigner function as 
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W(x,p;t) 



1 



2hir 2 
1 

4hn 5 



d 2 X 



; -A»7*+A*t) y r 



Pm (t) e 



Ab f -A*6 



d 2 (/3,/3,/3',/3')e- 2 '"'"* 



ie -PX-h\M 2 -k\P'X+Wc-n'+Vi' € '~ri +a2e -0i0'"-h\Pi?-h\P'l\ 2 +Wv*+Wi'n 



(V.2) 



Explicit expressions for ai,a2,«3 and a^, together with 
calculational details, are reported in Appendix D. 

Results for k = 2 and uj m — w™ p x 10~ 5 and 
cj m = 5w™ p x 10~ 7 , are shown in Figures [5] and [3] re- 
spectively for the homogeneous case. As seen for visibil- 
ity, in the case of (proposed) experimental values, Wigner 



function after measurement is undistinguishable from the 
free case. For the second (presently unattainable) much 
smaller value of w m , after a certain time a diminution of 
interference fringes is observed together with a contex- 
tual lowering of the first rest peak. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper the output of the mirror experiment pro- 
posed by Marshall et al. [29[ has been calculated within 
the framework of NNG, assuming both homogeneous and 
granular mass distributions. By varying the experimental 
parameters in a wide range beyond the proposed values 
[29| . a window of "sensible" parameters has been found 
in which the NNG induced decoherence effect is manifest. 

In conclusion, even if the experimental test of nonuni- 
tary gravity has been proved to lie beyond current tech- 
nology yet, requiring an unprecedented control of deco- 
herence, its peculiar form of self-gravitational interaction 
has been shown to be in principle distinguishable from 
the action of other collapse models. An exploration of 
the relevant parameter space could in fact, in a feasible 
experiment, lead to a clear distinction of the most appro- 
priate model. The signature of NNG model is ultimately 
connected with the fact that fundamental interaction oc- 
curs with a "simple" system (the mirror's copy in this 
case) rather than with a fundamental "thermal bath" 
random field leading to a visibility output somehow in- 
distinguishable from the effect of temperature. This is 
essentially due to the laboratory artificially created su- 
perposition state, while in naturally occurring circum- 
stances it is expected that "fundamental environment" , 
being as complex as the system itself, could eventually 
lead to auto-thermalization effects. 

In spite of the huge technical challenges, however, we be- 
lieve that due to the rapid progress in developing high- 
quality micro-optomechanical devices, a prototypal ex- 
periment of this type could be soon realized. 

A remark is finally in order concerning a finite temper- 
ature inclusion into the model. It should be clear that, 
when our initial knowledge of the system state is charac- 
terized by a density matrix like a thermal state, there is 
no unique prescription to associate it with a pure meta- 
state. In such a case one has to consider the possibility of 



using mixed meta-states to encode our incomplete knowl- 
edge. This more general case, independently of the spe- 
cific experiment treated here, will pave the way towards 
a generalized model of gravity induced thermalization. 
Such an interesting issue will be addressed in a future 
publication (37j | . 
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Appendix A: Interacting gravitational potential of 
the meta-mirrors 

In this Appendix the gravitational interaction poten- 
tial between meta-mirrors is computed within the two 
assumption of homogeneous and granular mass distribu- 
tions. 



1. The case with homogeneous masses 

Let's consider the Newtonian potential energy for two 
particles with masses Mi and M 2 



Starting from the initial condition where the two mir- 
rors of length L are overlapped, we consider the shifting 
of the first mirror of d/2 along Z\ positive axis and the 
second one along Z2 negative axis. Because of the very 
small relative displacement d between the meta-mirrors 



(a) 



(a) 




(c) 




(d) 





FIG. 2: (Color online) Wigner function in the homogeneous 
case for k = 2, cu m = wfn P x 10 _ , for the mirror size 
L — 10 m and for different intermediate times in a complete 
mirror oscillation. The variables x and p are in the ranges 

{-10Sx, 105x} and {-W5p, WSp} with 8x = 



(b) 




FIG. 3: (Color online) Wigner function for u„ 



hut 



10" 



All the other parameters are as above. All quantities 



and 



are expressed in dimensionless units. 



8p = 



while r = u*t. In the case with uj m = col 



as in the absence of gravity, we get qualitatively the same 
graphics. All quantities are expressed in dimensionless units. 
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(which is at maximum of the order of the size of the quadratic term of the expansion in the distance of the 
wave-function describing, in the ordinary setting, the CM total gravitational interaction energy 
coordinate of the mirror), it is enough to calculate the 



V(d) 



■I. 

-a 



(L-d)/2 AL+d)/2 

dz\ I dz 2 V (zi,z 2 ) 

J-[ 



(L+d)/2 
L/2 



-{L-d)/2 



dz\ + 



L/2 



Const. + - 




(L-d)/2 



dz 2 



(L+d)/2 



dz 2 V (zi,z 2 ) 



L/2 



where V (zi, z 2 ) 5z\5z 2 is the interaction energy of two 
square infinitesimally tiny mirror sheets parallel to the 
x — y plane, 



V(z u z 2 ) 



dxidyi / / dx 2 dy 2 - 



r\2 



Here we have used the assumption that, for example, 

I-( L L+d)/2 dz l S-L/2 dz ^(zi, Z 2 ) = f J^[ 2 /2 V(-%,Z 2 )dz 2 . 

In this way the terms linear in d vanish. Defining the 
non-dimensional coordinates x\ = X\/L, y[ = yi/L 7 and 
so on, we get: 



, L L\ ( L L 

vl r-2 )- v 2' 2 



G M 2 



1/2 



1/2 



dx 1 dy l I I dx 2 dy 2 



GM 
L 3 



-1/2 
2 A 



o J o 



-1/2 

i-e- 



1/2 



(A - 4) 2 + (y'i y'2) 2 + 1 



1/2 



d£ + / dr\. 
Jo Jo 



l-rj_ 



dr/+ 



1 



[e+v 2 -] 1/2 [£+v*_ + i\ 

1 1 \ 2irGM 2 



1 1/2 



G M 2 f 1 f 1 

/ / rf/?. (l-e_)(l-/7-) 



L 3 



Jo 



[e+v 2 



,1/2 



[e + v 2 - + 1 



1/2 



3 L 3 ' 



where we have introduced the new variables £± = x\ ± x' 2 gravitational potential between the meta-crystals de- 



and ?7 ± = y 1 ± y 2 . 
Then we obtain 



scribed by the state (to be symmetrized) 



V (d) = Const. + -GMp sil d 2 , 
o 

from which, writing the interaction term in the meta- 
Hamiltonian as — Kq (6 + 6^) (b + W), we get 



irhGp s 
3w m 



2. The case with matter-granularity effect 

Looking at the meta-mirrors as aggregates of atoms 
disposed in a lattice {Rj}, we are led to consider the 



i-d/2) \ 
Crystal J 



£(+d/2) \ 
Crystal / 



where superscripts parameters Td/2 refer to the center of 
mass displacements from the origin of coordinate system 
along, say, x— axis. The problem is similar to that of 
potential between two atoms or molecules, where (within 
the Born-Oppenheimer approximation) nuclei positions 
are treated as parameters of the atoms/molecules. 
The interaction potential is then given by (d = (d, 0, 0)) 
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V (d) = -Gml 



-Gmi 



-Gmi 



II 



dxdy 



y (-d/2) 
Crystal 



(x)V(x) 



(-d/2) 



Crystal I \ Crystal 



(+d/2) 



<A f (y)^(y) 



-,(+d/2) 

" Crystal 



( 



^nuc^Crystal \ 



N„ 



Hit 



E 

h,k=l 



dxdy- 



( m nuc LUCrystal\ , T f f , , < 

( M ) Nnuc J J dxdy - 



l x -y| 

i^-iJ( x _ R (-/ 2 )) 2 + ( y _ R ( + -/^)) 2 ] 

^£t£l[ (x+d /2) 2 + (y-d/2) 2 ] 



where, for simplicity, Einstein model for the crystal has 
been used, with uicrystai = Einstein frequency ~ 10 THz, 
m nuc ~ 4.7 x 10~ 26 Kg is nucleus mass, N nuc ~ 10 14 the 
number of nuclei. In the last line of the above formula, 
we have made the assumption that the nuclei wavefunc- 
tion spreads are much lower than interatomic distance, 
and that the formers are greater than d. To get a simple 
estimate, we can consider the interaction between inter- 
penetrating spheres of radius 2m '- ■ 
by a distance d, so that 



separated 



V (d) ~ Const. + -GMp nuc d 2 



(A.l) 



giving the (enhanced) gravitational coupling 



Kf an = 2HGpnuc ~ io 4 x K G om . 



Appendix B: Dynamical evolution of the 
meta-system 

In the following we study in detail the dynamical evo- 
lution of the meta-system. As initial meta-state we take 



I* (0)» = 2 KIOaIb) + |U0 B )) \a)] ® [(\0 A 1 S ) + 1^0^)) |S>] = |V> (0)) ® \j> (0)) 



(B.l) 



r 



where general coherent states a\ and ct2 are considered, Let's start by calculating the meta-state at time t. We 

although acceptable meta-states must be symmetrized will use the interaction picture, defining 
with respect to the physical and hidden parts. 

I 



| vj> (t)))j t = e i ( H /ree[6,6^X4,AAB;w^]+H /ree [6,6 t ,AA2',JV'B;w^])t/fi ||\I>(i))) 

I 



Sch 



where 



and 



\\*(t)))sch = \i>(t))®W(t)} 



j 



|^(t)) = ^e---*(|0^1 s ) 



ae 



g iK 2 (w^ l t-sina;^t) e iK9[a(l-e" i ^™ t )] |]^0 B ^ 
I 



ae 



C* + K (l- e~ iu ^ 



mi 



The interaction Hamiltonian 

H G = -K G (b + tf)(b + tf) 



evolves as 
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H GJnt (t) = _X Ge <^ee[6,fc t ^,AA B ;<]+ff /ra [6,6t^^~;<]) t /^ 5 + 6 t^ x 
X (6 + 6t) e -»( ff /rec[fc,b t ,A/'A,AAB;w^] + ff /ree [fci t ,A/'^,AA B J ;w^])t/?» _ 

= -K G e lHfr "\- b - bi ' -MA-NB^' m ]t/h ^ + ^ e -iH lr „[b,b\N A MB-M* m ]t/h x 

X e «H/ree[fc,6 t ,A/'A , A^b ; ] t /ft ^ _|_ >| g -iff /re e [i>, fc f , ^ , A/'g ;W^]t/ft_ 

I 



The above expression is the product of two specular and 
terms, so it suffices to calculate the first, say. We make 
use of the Backer-Hausdorff lemma: 

oo /_,Nn 

e- F Ge F = £ [ —r [F, G] n , F = -i [u,* m tfb gN A (b + 6+)] t, G = b + b\ 

n=0 n - 

with 

[F, G] = G, [F, G] n = [F, [F, G] n _ 1 ] , (g is now defined as g — klu^) by which, noting that 
I 

[F,G} = G = b + b\ 

[F, G], = -iLu* m t [6+6, 6 + 6+] = icj* m t (b - 6+) , 

[F,G] 2 = [F,[F,G] 1 ] = -{u* m tf[{b + b^)-2{g/u* m )M A ], 

[F, G} 3 = -*«t) 3 (6 - 6t) = -{u* m tf [F, G] 1 , 

[F,G] 4 = -(^t) 2 [F,G] 2 



\FA mm¥ a = (-l) n - n/2 (<t) n -" [(b + b^) - 2(g/u;*JM A \ 
[F,G] nodd = (^r-(6-6t), 

we obtain: 

O (t) = e iHfr " [ a > at >- /V 'A,A/'B;"i^]t/ft ^ -)- fttN g- i - H /rec[a,a t ,A/'A,AAB;w^]t/ft 

= 6 T (cos(w m i) + z sm(w m t)) + 6(cos(w m i) - i sm(u m t)) + 



OA™ UJ„ 



and a similar result for 0(t) : 



]t/ft 



= 6 T (cos(w m t) +tsm(u m t)) + 6(cos(u; m i) -%sm.{w m t)) — + — — , 

or 

6(t) = 7 (t)6t + 7 *( t ) 6 + r(t), d(t) = 7 (# + 7 *(t)6 + f (t), 

where 

7 (i) = cos(w^i) + isin(w^t), = -2/<jV J 4(cos(u;^i) - 1), f (i) = -2«A/A(cos(w^t) - 1). 

The evolution equation is 

mi ^ ))lnt = - l -H GJnt (t)\\*(t))) Int = l -K G 0(t)0(t)\\V{t))) InU 
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i.e., using the Hubbard- Stratonovich transformation [38j |. 



\\V(t))) Int = Teti kK G d{t')0{t')dt' ||^( ^ = 
x ( J 0|Te ic ^/ n -/'o d *'^ 1 -* ,2 ) 5 W|5)| ( 8>|ent} |3)|ent) = \ If d 2 pd 2 pK UA * A {(3,(5- a -t)\(3)\ent) ® |0}|ent>, 



where c is a constant. Before dealing with the kernel 
j^MaMa ^ we calculate the amplitude with the help of the 
Bakcr-Campbell-Hausdorff formula: 

e t(A+S) = e tA e tB e -4[A t B] e ^(2[B,[A,B]]+lA,[AM) 



= (p\Texp(fbt f di'(^i + ^ 2 )7(«)}exp|/6y + </? 2 ) 7 *(i)} x 

x exp{^(7 dffai + pafrwV^ ^(^+^ 2 ) 7 *(t))}|a)exp|/y df fai + pa)r(t)} = 

= e /^* /o dt'( ¥ >i+ V2 )7(t) e /a/ t rft'( V i+¥> 2 )7*(t) e ^ (/o <«'(¥>i+y2)7(t)) (/o dt'tei+^h'W) g// * dt'( yi + V2 )r(t) 



A similar result holds on for the other am- 
plitude by making the following substitutions: 



We get then 
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~ r r r err ( ' 2 JV 

K N A N A (j9j ^ a; t ) = Jim //.../ JJ #2ii ... ft d Vl , exp -^f £ A *K* - 

At^O JJ J »=1 ^ »=1 i=l 

JV 2 JV . 

+ /£Ai(<^ M + ^ 2 , i )(/r 7i + a 7 *+r i ) + — 2(At)Vi,i+^,ihi7i(Pij+^j)f x 

i=l i,j=l ^ 

x exp<^ /^At( Vlii - <P2,i){p*ii + a 7 * + y ( A *) 2 (^M - V2,i)7i7;(Vi,j - ^2,,) [(/3|a)(/3|5) = 

i=l i,i'=l ^ 



= exp 



x lira 

JV^co 
Ai^O 



lal 2 \B\ 2 lal 2 ~ 

^L+/3*a-^--^-+/3*a 



, JV 



JV 



^A^ J + ^(At) 2 7j7 ; 



//•••/ n^^^pjn 

2—1 i,i= 1 



7C 



/ 2 



_ A ^ + T -(A^ 7i7 ; 



V2,^2,j + /£ A*(/r 7i + r * - z 3 *^ - r 4 )^ 2jl ] X 
i=i > 

JV 

<Pi,i<Pi,j + /£ A*(/r 7i + 2a 7 * + r, + /?* 7j + r\> v , 



= e 2 2 



+/3*" e -^p-^+/3*« x i im 



2 i „. 1 2 ~ 



(27T) 



2JV 



jv^oo V det AW 

At^O 



x lim 



(27T) 



e 2 



JV 



ij^AW)"^ = 



te^/dctAWdetAW 



ijr(A< N) )- i Ji-ijr(A< N ')- i J 2 



A(t) 



x e -§l /o/o ^"[/^(t^a^tO+rcO+zS^ x 



x e 2R 



^ /J/ * dt'dt'qr 7 (0+r(tV/3*7(0-r(0]^^^^ 



where 



¥>1,1 = <Pl(t = 0), ¥7l,JV = ^2,1 = ^2^ = 0), l^ 2 ,JV = <f2(i); At=—, T = U* m t, 



and 



A (Jv )= /AW 0. 



A 



(JV) 



j(JV) = ( Jl 



(JV) 



,*-2 



(At)' 



ft At ft 2 -"IT. 7jj 



AN) 



■ *— 2 



(At)' 



ft At ft 2 ^LT. 7,j 



(JV) 



(r +/?*) 74 + 2a 7 *+r^ 



(/3*-/3*) 7i + r^ 



here with 5ft [x] we indicate the real part of x and We take now the continuous limit as 

t m a kt a = r ( t ) + f ( t ) and r^- 4 ^ = r(t) - f (t). 
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lo* 2 A (N) 



J<"> 



At->0 



At at^co 

•J_2 

At at^oo 

At->0 



A 1[2] (t',t" 

Ji(r) = 
J 2 (r) = 



1 



K G 
Kg~ 



A = lim 



(27T) 



AT 



jy^oo c 2 vdetAi dct A 2 

Ar->0 



with c 2 



At 



Alternatively, the function A(t) can also be obtained from 
the normalization condition ((^(t)\ = 1. 
Defining now, for general functions /(t) and g(r), 

T T 

4?]to = / ldT'dr"f(T')A- [ l ] (T',T")g(T"), (B.3) 





08* + /3*) 7 (t) + 2a 7 *(T) + P^(t) 

(^*-,s*)7(T) + r^^ A (T) 



r 



m we get 



A(t) e -T-^+P'" e -T^ 



|/3| 2 _|a| 2 ~- 



X 



l^C "f [(r-r ) 7 (-')+ r - A ^ (-')] A ^ (t'.t") [(r-r )7(r")+r^^ (r")] 



x e 



= A(i) e 2 « e 2 —+P « 



x e 



x e 



x e 



x e 



x e 



x e 
Finally 



'+fS') 2 f f dr'dr" 7 (r')A 1 - 1 (T',r") 7 (r")+4a(/3* + /3*) / / dr'dr" 7 * (r') A" 1 (r',r") 7 (r") 



2(/3*+/3*) / / dr'dr" 7 (- 



-') A^Cr'.r'Or^Cr")! 



3*-/3*) 2 "/ ""j t dr'<ir" 7 (T')A- 1 (r',T") 7 (r")+2(/3--/3*) "j^fdr'HT 1 ^^') AJ^t'.t") r^^(r") 



4a7*TdT'*"'V(T') ^OV") <^(r")J 



X 

t UJ* t 



4a 2 "/ *"/ t dr'dT" 7 *(T')A 1 - 1 (T',T") 7 *(r")+ / / dr'dr'T^^^ (t') A" 1 (r',r")r 





^(r»)} 



/ / dr'dr'T^^^^^A- 1 (r',r") T^ A " A (r") 



K 



•V,V„ , -, n .. V „ , „-J^-J^ + ^ Q+ ?. a -| Q | % -^^{^[^ ,r+ ^^*] + ^ + ' r+ + ^- r - } 



/?; = A (i) e 2 2— «-r^ e i. L - j - Jx 

i T /f G {^7^(^*+r) 2 +^^(^*-r) 2 +2(/3*+r)[2 Q ^7*^ +y;- r +]+2^- r -(p*~0') } 

I 



(B.4) 



x e 



Inverse operators Aj^ 1 , A 2 1 have been evaluated numer- Note that the dependence on /?, /3 is present only in the 



ically. 



second factor (second and third row), while the depen- 
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dence on photon states is hidden in T^ aMa , T^ aMa . Appendix C: Integrals of visibility 

Schrodinger state at time t is then given by Equation 
PV.1| ) • Integrals (I) and (II) appearing in Equation (|IV.2[) arc 

given by 



(J) = _L e « 2 K*„t-si„ <t) K io (i) K oo* ^ x 



(II) = J-e^^-^^K 11 (t)K al * (t) x 

4-7T 4 

x || dWe"^ 1 "^^^ (/?,/?') h((3, (S')Kl 1 {p) Kf (13') J f <?p<Fp> ■K^ 1 0)I^ 1 ^) e 0)S(^)g0,?)H l (fi,p'). 

I 

Appendix D: Calculation of the Wigner Function p m (t) . Defining 

In the following we present the details of the calcula- 
tion of the Wigner function. Let's start by computing 



II* (*)»«, = M®MII*(t)» 



1 



d 2 (3d 2 (3 



if uu (/3,/3) |/3 C ) /3 C ) + K^(p,(3)f(P)e- l » \p c ) 



A) + 



I 



and 



the density matrix of the physical mirror is: 



pS) = \m))um)\\, 



Pm {t) 



Tr(p v (t)) = - / d'x (X\ P V \X) 

7T 



4tt 4 



d 2 (/3, /3, /?', /?') ai |&) (#| + a 2 |A> <#| + «3 |&) (A'l + "4 \Pi) iX\ , 



where 



ai(j3j,/3'J') = K 00 (p,p)K* 00 (p\p')e^''-^\ 2 -^\ 2 + 

+ K 00 (f3j)K*°\(3 / J')e td f%^)e^ J ''-^\ 2 -^\ 2 + 

+ K ol (^l3)K* o \p\P')f0)f*0')e^''-^'\ 2 -^''\ 2 , 
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a 2 (M,p'J') = K l0 {j3,P)K* 10 {p',p')f{P)f*{P') e ^'""-i\^\ 2 -^\ 2 + 

+ iif 10 (/3, ^) J ftT* 11 (/3',/3')/(/3)/*(/3')e ie /*(/3') e^^'-i l^l 2 -!!^^ I" + 
+ K 1 \(3j)K* lo (P\P')e- ie f(P)f^(3')f0)e^'"~^ l \ 2 -^\ 2 + 



+ ^ oo (/3J)K* 11 (/3^/30e 2^ V*(/30^(?')e^*"^^| 2 "^^| 2 + 

+ K o \pJ)K* 1Q {f3',f3 / )f0)f*((3')e^''-^ l \ 2 -^\ 2 + 

+ X ol (^J)if* 11 (/3^^0e ^ V(/3)/*(/30^(/?')e^^*"5l*| 2 "5l^l^ 



a 4 (p,j3,p'J') = K w ([3j)K* 00 (P\p / )e- ie f(P)e^~^\ 2 ~i\P'X + 

+ K lo (l3,j3)K* ol (P'J')f(P)f*0')e^-^\ 2 -^'^ + 

+ K n ^J)K* oo (/3\^)e- 2i6 f(/3)f0)e^'*-i\M 2 -M\ 2 + 

+ K 1 \pJ)K* ol ^'J')e- ie f(p)f(j3)f*0')e^^-^ l \ 2 -^\ 



By taking the trace 



Tr 



Pm (t) e 



- I d 2 X { X \p m (t)e xb '- x ' b \x) 



we finally obtain Equation (|V.2[) . 
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